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ABSTRACT 


Equations are derived for maximizing payload capability for nonplanar boost trajec- 
tories arising from range safety considerations. The assumed range safety constraints 
are satisfied by requiring a fixed vehicle azimuth heading during the initial part of the 
trajectory. This azimuth (attitude) constraint is removed as soon as certain land areas 
have been avoided. An attitude rate limit is also imposed to prevent discontinuous atti- 
tude changes. The calculus of variations is used to formulate and solve the mathematical 
problem. A simplified method is also presented for obtaining the auxiliary variational 
terminal conditions associated with free end conditions. Sample results are presented 
for sun-synchronous and planetary missions . 
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SUMMARY 

Equations are derived for maximizing payload capability for nonplanar boost trajec- 
tories arising from range safety considerations. The assumed range safety constraints 
are satisfied by requiring a fixed vehicle azimuth heading during the initial part of the 
trajectory. This azimuth (attitude) constraint is removed as soon as certain land areas 
have been avoided. An attitude rate limit is also imposed to prevent discontinuous 
attitude changes. The calculus of variations is used to formulate and solve the mathe- 
matical problem. A simplified method is also presented for obtaining the auxiliary vari- 
ational terminal conditions associated with free end conditions . Sample results are pre- 
sented for sun-synchronous and planetary missions. 


INTRODUCTION 

Injection conditions for most earth orbital and planetary missions are achieved by 
using essentially planar steering throughout the entire launch vehicle-powered flight tra- 
jectory. A planar trajectory results from selecting the launch azimuth heading which 
gives the desired orbital inclination or declination and right ascension of the outgoing 
asymptote without yaw steering. However, range safety instantaneous impact point (IIP) 
restrictions sometimes do not permit the use of the launch azimuth required for a planar 
trajectory. In such cases, the required injection conditions may be achieved by employ- 
ing a nonplanar trajectory and steering the vehicle so that the probablility of impacting 
land masses is eliminated or at least minimized. 

The rigorous optimization of trajectories with HP constraints requires a detailed 
specification of the IIP restrictions. These restrictions must then be formulated mathe- 
matically and included along with the vehicle equations of motion and mission constraints 
and requirements. This problem can be solved by using the calculus of variations. How- 


ever, the IIP restrictions result in an extremely tedious mathematical formulation. 
Therefore, a simplified solution procedure is presented in this report. A launch azimuth 
heading is chosen which allows the initial portion of the trajectory to avoid all land 
masses. The thrust direction is restricted to be parallel to the launch azimuth plane (the 
direction within the plane is optimized) until the HP has passed the protected land masses. 
Then, the planar steering constraint is removed and optimal three-dimensional steering 
is used to satisfy the required injection conditions at powered flight termination. 

At the switching point between planar and unrestricted steering, the optimum trajec- 
tory requires a discontinuous change in vehicle attitude, if no attitude rate limits are 
imposed. Since a discontinous attitude change is not physically possible, an attitude rate 
limit is added to the problem formulation. 

The calculus of variations has been used extensively in the past to solve both two- 
and three-dimensional trajectory problems. Reference 1 derives equations for optimal 
two-dimensional ascent trajectories and also includes equations for optimizing stage pro- 
pellant loadings. Equations for optimal three-dimensional interplanetary trajectories are 
derived in reference 2. Equations have been developed by Johnson and Winfield (ref. 3) 
for rate limited optimal trajectories. In that report, the rate limit problem arose because 
of constraints on initial and/or final thrust direction. Consequently, the rate -limited 
phases (referred to as regular arcs) occur at the beginning and/or end of the trajectory. 
Intermediate attitude rates are used throughout the rest of the trajectory. These inter- 
mediate attitude rates result from the existence of singular arcs in the optimal trajectory. 
The equations are derived by using the maximum principle, and a possible solution pro- 
cedure is discussed. However, no numerical examples are solved. 

The present report extends the analysis of reference 3 to include constraints on 
vehicle attitude, as well as attitude rate. Also, the problem solution requires regular 
arcs in the center, rather than at the ends, of the trajectory. The calculus of variations 
is used to obtain the mathematical problem formulation. A simplified method is also 
presented for deriving the auxiliary variational boundary conditions associated with free 
trajectory injection conditions. 

In order to illustrate the use of the equations, solutions are presented for two typical 
missions which require nonplanar trajectories of the type described. An Atlas-Centaur 
launch vehicle is assumed for both of these missions with an Eastern Test Range (ETR) 
launch. The first mission is a high altitude (1110 km) sun- synchronous orbit, with a 
required orbital inclination of 99. 9°. The second mission presented is a typical Mars or 
Venus probe with a declination of the outgoing asymptote of -45°. 


ANALYSIS 


The problem to be solved is to determine the optimal thrust direction profile for a 
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trajectory satisfying certain boundary conditions. The thrust direction is constrained to 
lie in a fixed inertial plane during the initial part of the trajectory; later, the thrust direc- 
tion is completely free in three dimensions. However, the attitude rate is constrained to 
be less than or equal to a prespecified limit throughout the variational portion of the tra- 
jectory. At the transition point between the constrained and unconstrained attitude por- 
tions, a discontinuity in thrust direction occurs in the optimum trajectory if no attitude 
rate limit is imposed. However, a maximum turning rate is specified in the present 
problem so that the vehicle will rotate at the maximum rate (regular arc) in this region. 
The boundary value discussion presented in this report does not consider the possibility 
of regular arcs occurring in other regions of the trajectory. Generally, the maximum 
turning rate is large enough so that additional regular arcs are not required. If additional 
regular arcs do occur, the boundary value problem must be altered accordingly. 

Problems of this type can be formulated as a Bolza problem (ref. 4). Specifically, 
the functional to be minimized is written as : 


J = 



(la) 


where 


and 


F = X • bj + ii ■ b 2 + cr • bg + yb 4 + pbg 


(lb) 


g = t f (lc) 

subject to the following constraints: 


L = v+ — r - af = 0 
1 3 

(2a) 
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(2b) 

bg=f-COXf = 0 

(2c) 

b 4 = w • w - = 0 

(2d) 

bg = (3J x f) • n = 0 

(2e) 
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All symbols used herein are defined in appendix A. 

Equations (2a) and (2b) are the vehicle equations of motion written in three- 
dimensional rectangular coordinates. The thrust acceleration profile, a, is assumed to 
be a definite function of time . The unit vector f defines the thrust direction. Equa- 
tion (2c) determines the rotation rate I as a rotation of f about the vector u>. Since 
f is perpendicular to f , the magnitude of f is constant (as required) and is set equal to 
unity as an initial condition. The optimum vehicle turning vector, a>, will be determined 
by the variational analysis. However, the magnitude of f must not exceed the maximum 
turning rate, a> m . This constraint is satisfied by equation (2d), which states that the 
magnitude of o> is equal to w m . From equation (2c), the magnitude of t is equal to 
w m ^ k* is perpendicular to f; otherwise, the rotation rate is less than a> m . Since no 
restriction has been placed on the direction of aJ, equations (2c) and (2d) are equivalent 

A 

to an inequality constraint on f . 

Equation (2e) is used only when the planar attitude constraint is in effect. (The unit 
vector n is the normal to the required thrusting plane.) When the attitude is uncon- 
strained, p is set equal to zero. Equation (2e) requires the component of i along n to 
be zero (i.e., f • n is constant). The value of f • n is set equal to zero as an initial 
condition. 

The variables X, p, a, y, and p are undetermined Lagrange multipliers which are 
functions of time since the constraint equations must be satisfied at all points of the tra- 
jectory. The Euler- Lagrange equations for this problem (ref. 4) are 


X + p = 0 

(3a) 

£-£H- + 3S2(a. 7)7=5 

(3b) 

3 5 


r r 


q + aX+qXa> + pw x n = T5 

(3c) 

2ya> - fxq+pfxn=TJ 

(3d) 


The turning rates for the planar and unconstrained attitude cases must be considered 
separately. For the unconstrained attitude case (p = 0), combining equations (3d) and (2d) 
results in 


co 


= ±co 


m 


? X q 
f X q 


(4) 


Thus, the optimal turning rate has been determined, except for sign. The proper choice 
of sign will be determined later. Since w • ? = 0 (from eq. (4)), the maximum turning 
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rate (regular arc) is always used, unless f x a =TJ in which case co is indeterminate. 
It will be shown later that f X a = 0 corresponds to the singular arc. 

For the planar attitude case, equations (2d) and (3d) give: 


CO = ±C0 

m 


f x g - pf x n 

f x a - pf X n | 


The value of p is determined by combining equations (2e) and (3d). This procedure 
results in: 


p = (f x a) • (f x n) 

If f x cr is expressed in the orthonormal f, n, f x n coordinate system, 
f x a = [(I x a) • f]f + [(f x a) • njn. + [(f x a) • (f x n)jl x n 
and since (f Xcr) • f = 0, 

lx? - [(f x o) • (f x n)]f x n = [(f x a) ■ n]n 


or 


f x g - pf x n = [(i x a) • njn 


Therefore, 


co = ±co m sign [« x a) • fi]n (5) 

Since a> • f = 0, the maximum turning rate is used on the regular planar arc, just as on 
the regular unconstrained arc. The singular planar arc occurs when (J x <j) • n = 0. 


Weierstrass Condition 

The uncertainty in sign in equations (4) and (5) can be resolved with the aid of the 
Weierstrass E-test (ref. 5), which can be stated as 

E ^ 0 for a minimum 
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where 


E = F(x* xp - F(x i , 



( 6 ) 


The x^ correspond to the minimizing values of the problem variables (r, v, f, and To), 
which differ from the x* by a finite but admissible amount. For the present problem, 
only To is subject to such a variation, and equation (6) can be written explicitly as 


E = a ‘ (co - co*) x f > 0 for a minimum. 


Since o>* = -co, this equation becomes 

To ■ (a X f) < 0 


(7) 


for a minimum. 

For the unconstrained attitude case, combining equations (4) and (7) results in 


±w. 


m 


(? X q) • (q X j) 

If X q| 


< 0 


or 


±<X)„ f X q >0 
m 1 


which requires that the plus sign be chosen in equation (4): 

f x q 


co = co 


m 


f X q 


For the planar attitude case, equations (5) and (7) are combined to give: 

±o> m signj(f x q) • njn • (?xf) <0 


or 

±w m |(f x a) • n| 2= 0 

Again, the plus sign must be chosen. Therefore, 


( 8 ) 
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(9) 


0 ) 


= sign 


p X cr) • njn 


Singular Arcs 

Singular arcs occur when the Weierstrass E-function is identically zero for some 
finite time interval. From equations (4) and (7), this occurs for the unconstrained atti- 
tude case when f x a = 0. If f x cr = 0 for some finite time, then (d/dt)(? x a) must also 
be zero. But from equations (2c) and (3c) (with p = 0), 

— (fXa) = (aJxf)xa+fx Jj-aA - ct x cJj 
dt 


= (ojxf)Xcr-afxA-fx(aXco) 


= -af xl - wx(ffxf) = -af x A = 0 


( 10 ) 


which implies that 


f = ±- (11) 

A 

for the singular arc. This is the same solution as obtained in reference 2, where no rate 
limit is imposed. It is shown in reference 2 that the Weierstrass E-test requires that 
the plus sign be chosen in equation (11). Therefore, 

f=A (12) 

A 

is the singular solution, which takes effect when equation (12) and f x a = 0 are satisfied 
simultaneously. 

For the planar attitude case , the singular arc and corresponding indeterminacy in 
uj occur when (f X a) • n = 0 (eq. (9)). The required thrust direction for the singular 
arc is obtained from 
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( 13 ) 


— (f x a) • n = (a) x f) x a • n+fx £-aX - a x a> - pa> x n] • n 
dt 

= (w x f) x a • n - af x X • n - f x (g x w) • n-pfx(wxn) • n 
= -a(f XX)- n = 0 

Since f is normal to n and X X n, f may be written 

t = n x (x x n) 

| x x n| 


or 

f _ X - (X • n)n (j 

Vx 2 - (X - S) 2 

where the sign choice in equation (14) has again been made by using the Weierstrass E- 
test on the unrestricted rate, planar attitude solution. 


Weierstrass-Erdmann Corner Condition 

The boundary conditions on the Lagrange multipliers at corner points of the trajec- 
tory can be derived with the aid of the Weierstrass-Erdmann corner condition (ref. 4). 
This condition requires that the 3F/ Sx^ be continuous at all corner points of the trajec- 
tory, where the are the problem variables as defined earlier. For this problem, 
corners occur at transition points between regular and singular arcs and at the switchover 
point between planar and unrestricted attitude. The 3F/3x^ are X, p, and a; therefore, 
these variables must be continuous throughout the entire trajectory. 


Trajectory Description 

A list of major trajectory event times is presented in table I, along with a brief 
description of each event. The trajectory starts with a short vertical rise, followed by 
a rapid pitchover phase in the desired azimuth direction, which determines the amount 
of trajectory lofting during the atmospheric phase of the trajectory, as well as the vehicle 
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azimuth heading. The rest of the atmospheric phase (which is assumed to be terminated 
at booster stage burnout) is flown with a near-zero angle of attack steering program (such 
as described in ref. 6), in order to minimize vehicle heating and aerodynamic loads. The 
thrust direction is maintained perpendicular to n throughout the booster portion of the 
trajectory. 

The variational problem begins after the sensible atmosphere has been traversed. It 
is assumed herein that this occurs at fixed time, position, and velocity. This implies 
that the magnitude of the initial pitchover is fixed. However, the degree of freedom avail- 
able when the pitchover magnitude is free for optimization can also be handled by varia- 
tional methods, as shown in references 1 and 6. It is also assumed that the desired 
initial thrust direction for the variational phase can be achieved prior to the start of this 
phase, so that an initial regular arc is not required. This assumption is reasonable since 
the attitude discontinuity at this point is generally small (see Results). The equation 
modifications required to handle an initial regular arc are straightforward and are 
described in reference 3. The thrust direction during the initial portion of the variational 
trajectory is constrained to be normal to n, as in the booster trajectory. Thus, a 
singular, planar arc is used. 

The initial singular, planar arc continues until a short time prior to the time at 
which the planar steering constraint is removed. At this point, a regular planar arc is 
used. This arc is not forced by saturation (i. e. , the variational commanded planar turn- 
ing rate does not exceed u> m ). However, a regular unconstrained arc is forced at the 
junction point between the planar and unconstrained portions of the trajectory, because 
the conditions for singularity are no longer satisfied at this point. The regular planar 
arc is therefore used in order to shorten the forthcoming regular unconstrained arc. In 
addition, the use of the regular planar arc supplies an additional degree of freedom (start 
time of regular planar arc) which will be needed to satisfy the boundary conditions to be 
discussed later. A second degree of freedom results from the fact that f x o is not com- 
pletely specified on the planar arc. Johnson and Winfield (ref. 3) discuss the fact that no 
additional degrees of freedom are available for regular arcs forced by saturation. How- 
ever, two additional boundary conditions must be satisfied for each regular arc intro- 
duced, and these boundary conditions can only be satisfied by anticipating the forthcoming 
regular arc, thus allowing two additional degrees of freedom. For this reason, regular 
arcs forced by saturation only exist as limiting solutions. 

The need for introducing a regular arc prior to saturation may be better understood 
by referring to a two-dimensional analogy, as illustrated in sketch a. 
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(a) 

Suppose the optimal pitch attitude history for a nonrate-limited trajectory is given by 
curve 1. If a rate limit is introduced, the trajectory may be modified by using a regular 
arc starting at the saturation point (curve 2) and terminating when the attitude coincides 
with curve 1. The performance loss associated with curve 2 is approximately propor- 
tional to the integral of the absolute attitude difference between curves 1 and 2. This 
integral can be decreased by initiating the regular arc earlier (curve 3) , thereby improv- 
ing performance. The use of curve 3 also introduces a degree of freedom (time of regu- 
lar arc initiation) , which is needed to satisfy a variational boundary condition introduced 
because of the regular arc. 

When the attitude constraint is removed, the trajectory switches from a regular 
planar arc to a regular unconstrained arc. This arc is then continued until the required 
boundary conditions for a singular arc are satisfied. At this point, a singular uncon- 
strained arc is used and is continued until the desired end conditions are achieved. It is 
assumed that the final attitude is unconstrained; therefore, no final regular arc is 
required. Again, the equation modifications required to include a final regular arc are 
straightforward and are discussed in reference 3. 

TRANSVERSALITY EQUATION 

The relation between changes in boundary values and changes in J (eq. (la)) is 
expressed by the general transversality equation (ref. 4). 
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where dJ has been set equal to zero for an optimum solution. For the present problem, 
dJ is evaluated to give 


d J = (C dt + X • dv + p • dr + o • df L + dtj, = 0 


where C is a constant, the first integral of the Euler-Lagrange equations (ref. 4). 


C = X - r-p-v-aX-f-cr - (c3 X t) 
3 


Since it has been assumed that the initial position, velocity, and time are specified, 
dr Q = dv Q = dt Q = 0 and equation (16) becomes 

[<C + 1) dt + X • dv + ju • dr + a ■ dj] tf - (5 • ctf) t = 0 (17) 

The final time is not specified so that C = -1 for an optimal solution. Since the Euler- 
Lagrange equations are homogeneous in the multipliers, C = -1 can be satisfied by the 
proper choice of multiplier scale factor. In any case, C does not appear in any of the 
equations used to optimize the trajectory. 

The final attitude is not specified, hence (a • df), = 0 must be satisfied for an opti- 

r f 

mal solution. The allowable variation in f (perpendicular to f since f is a unit vector) 
is given by 

A A 

df = f X 6 


where e is an arbitrary vector. Therefore, 


(a • df), = (a • f x e). = (a x f) • e. =0 


Since e is arbitrary, 


(ffXf), =0 
t f 


is a required boundary condition. 
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At t Q , f has only one degree of freedom since f • n must be zero. This degree of 
freedom is given by 


df = af x n 


where a is arbitrary. Combining this result with the last term in equation (17) gives 

(a x f), • n = 0 (19) 

l o 

for a required boundary condition. Equations (18) and (19) are also necessary conditions 
for singular unconstrained, and singular planar arcs, respectively. 


Boundary Conditions on Position and Velocity 

There exists an auxiliary variational boundary condition to be satisfied for each free 
trajectory end condition. By satisfying this auxiliary equation, the resulting value of the 
free end condition will minimize flight time . The auxiliary equations are obtained from 
the position and velocity terms in equation (17). 

X • d y + Jl • dr = 0 (20) 

Suppose the final orbit is specified in terms of arbitrary independent variables q^ 

<lk = < lk,d k = l,2,. . .,n^6 

If n is less than six, there are free end conditions. The set q^ is completed by addi- 
tional independent variables, if necessary, to give a complete specification of the tra- 
jectory state in terms of the q^. By expressing the components of r and v in terms 
of the q^, and differentiating, equation (20) becomes 

6 

Z — — (A • v + ju • r) dq^ = 0 

3q k 

k=l K 


For each free q^, the corresponding auxiliary boundary condition is obtained by setting 

— — (A • v + ju. • r) = 0 (21) 

3q k 
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This procedure, although straightforward, can be very tedious if the q's are com- 
plex functions of r and v. For this reason, a simplified approach is presented for cal- 
culating the required auxiliary final conditions . With this approach, the dr and dv are 
interpreted as allowable variations in r and v. Given a set of q^, the allowable vari- 
ations are such that all other boundary values remain constant, while the free end condi- 
tion is allowed to vary. In other words, the allowable variations correspond to partial 
derivatives of injection position and velocity with respect to the free end conditions, as in 
equation (21). To illustrate the use of this procedure, the auxiliary boundary equations 
will be derived for the two missions considered herein. 

Sun- synchronous orbit . - The required orbit conditions for the sun-synchronous orbit 
mission call for a circular orbit at a specified altitude (radius) with a specified orbital 
inclination. This is equivalent to the specification of four orbital parameters: 


q l = r = r d 



q^ = r = v • r = 0 


q 4 " 1 “ J d 

Two additional independent orbital parameters must be added to complete the set of q^, 
and these additional q^ are free for optimization. The free orbital parameters may be 
considered to be 


q^ = longitude of ascending node 
qg = argument of injection latitude 

For free ascending node, the trajectory plane may be rotated about a vector (z) 
pointing at the North Pole. Such a rotation changes the ascending node but leaves 
all other orbit parameters fixed. The allowable variations in position and velocity 
corersponding to this rotation are 


dr = ar x z 
dv = a v x z 
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where a is the arbitrary magnitude of the rotation. Combining these variations with 
equation (20) yields 


ot(\ • vXz + ju • rXz)=a(Ixv + jiXr)' z = 0 


Since a is arbitrary, 


(A x v + (i x r) • z = 0 


( 22 ) 


for optimum ascending node. 

For free argument of injection latitude, the allowable variations in r and v result 
from a rotation of the trajectory about the angular momentum vector, h. Hence, 

dr = ar x h 

dv = av x h 

and the auxiliary variational end condition is 

(XXv + juXr)- h = 0 (23) 

Interplanetary transfer. - The required heliocentric transfer trajectory may be spec- 
ified in terms of the required hyperbolic excess velocity vector, v M , at Earth injection. 
This vector, in turn, may be specified in terms of three variables: 

q^ = injection energy (or equivalently, hyperbolic velocity magnitude) 
qg = declination of outgoing asymptote 
q^ = right ascension of outgoing asymptote 

In addition, the perigee radius of the geocentric hyperbola is specified. 

«4 = r p 

As in the sun-synchronous mission, there are two free end conditions. These cor- 
respond to 

q^ = injection true anomaly, Vj 

qg = a free rotation of the trajectory plane about the hyperbolic 
excess velocity vector, v^ 
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For free true anomaly, the allowable dr and dv may be obtained as changes result- 
ing from coasting along the geocentric hyperbola at the injection point. Hence, 


dr = orv 


dr = qiG 


where 


G = - 


Gm- 


r 


3 


and a is arbitrary. Combining these variations with equation (20) results in 

X • G + p • v = 0 (24) 

for optimum true anomaly. 

The auxiliary variational end condition for qg is obtained in a manner similar to 
equations (22) and (23). This equation is 

(\ x v + n x r) • v w = 0 (25) 


NUMERICAL ITERATION 

An optimum trajectory requires the simultaneous integration of equations (2a), (2b), 
(2c), (3a), (3b), and (3c). As initial conditions, these equations require the initial posi- 
tion and velocity, as well as the multipliers X and /i. Since the variational trajectory 
begins on a singular arc, the value of f x cF is not needed at this point. The trajectory 
is then uniquely specified until the transition time to the regular planar arc (t^) . This 
transition time, as well as the component of f x a along f x n (f x a • n must be zero 
while on the singular planar arc), must be supplied as initial conditions at this point. 

The trajectory then proceeds along the regular planar arc and into the regular uncon- 
strained arc (at tg). On the converged trajectory, the transition to the unconstrained 
singular arc (tg) will occur when the equations for singularity are satisfied. During the 
iteration process, however, this transition occurs at a specified time (varied during the 
iteration process) and jump discontinuities in f and f x a to the singular values are 
permitted. The trajectory integration then proceeds uniquely to final injection. 
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With the initial conditions as given above, the trajectory is uniquely specified. It 
remains to find the values of the unknown initial conditions which result in the trajectory 
satisfying the required final conditions . The initial and final conditions for the two point 
boundary value problem are as follows: 


Unknown initial 
conditions: 


Desired final 
conditions : 


Mt 0 ) 

M(t 0 ) 

(f x a) • (lx n)(tj) 


Six orbital and/or auxiliary 
final conditions 

I x a = 0 (two conditions) 

£ • f = 1 

X 

C = -1 


(26) 


The boundary value problem may be simplified by using the initial value of one of the 
multipliers to satisfy C = -1. Also, the trajectory can be terminated on one of the final 
orbit conditions, thus eliminating tj and one final orbit condition from the problem. The 
resulting iteration size is eight by eight, and the iteration variables are 


Unknown initial 
conditions : 


Desired final 
conditions : 




J Mt 0 ) 

Five of < _ 

V<t 0 ) 

(f Xa) • (f xi)(tj) 


t 

t 


1 

3 


Five orbital and/or auxiliary 
final conditions 

I x a = 0 (two conditions) 

A . f = 1 

X 


(27) 


The iteration process is carried out by using a multivariable Newton- Raphson 
method. This method is described in detail in references 1 and 6 and will not be dis- 
cussed here. The only difference between the procedure described in these references 
and that used herein is the method of obtaining the required partial derivatives of final 
conditions with respect to initial conditions. In references 1 and 6 this was accomplished 
by using finite differences. This procedure was also attempted on the present problem. 
However, because of the sensitivity of the derivatives and the size of the iteration loop, 
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convergence was very difficult to attain. For this reason, it was decided to use integrated 
partial derivatives. With this procedure, convergence was easily obtained. The equations 
used in obtaining the required partial derivatives are presented in appendix B. A com- 
plete list of equations required to optimize the trajectory is also presented in this 
appendix. 


Choice of Initial Conditions 

In order to facilitate convergence of the two point boundary value problem, it is 
desirable to start the iteration with initial conditions which are as close as possible to 
their converged values. With this in mind, initial conditions have been chosen which have 
more intuitive significance than X and / u . The initial conditions chosen are presented 
in appendix C, along with the equations for calculating X and ju from these values. The 
rationale used for selecting initial guesses of the other unknown initial conditions is also 
presented in appendix C. 


RESULTS AND DISCUSSION 

The solution to the limited thrust attitude rate problem is most interesting in situ- 
ations where a significant period at maximum rate is required. One area where this 
occurs is when the instantaneous impact point (HP) constraints are factors in determining 
an acceptable trajectory. For launches from the Eastern Test Range (ETR), near polar 
orbits and hyperbolas with certain combinations of declination and right ascension of the 
outgoing asymptote require launch azimuths for planar trajectories which appear to be 
unacceptable because of the HP traces. Dog-legged trajectories may be used to try to 
improve the traces by launching at a nonoptimum azimuth and constraining the thrust 
attitude to be parallel to the launch azimuth plane until some later time in the trajectory. 
Two cases were chosen to demonstrate that solutions could be obtained with the analysis 
described even though the eight-by-eight iteration required appears formidable. The 
procedures used to alter the IIP traces in the cases are examples of techniques which 
might be used to improve the traces. However, the traces presented do not necessarily 
satisfy range safety constraints. 

The vehicle used in these demonstration cases is the Atlas -Centaur. The Atlas is 
propelled by two booster engines and one sustainer engine. The booster engines are 
jettisoned at a predetermined acceleration level. The sustainer engine continues to burn 
(sustainer solo). The Centaur insulation panels and then the payload fairing are jetti- 
soned in this phase. The sustainer solo ends at propellant depletion and the Atlas stage 


17 


is jettisoned. After approximately 10 seconds, the Centaur engines ignite and burn con- 
tinuously until the desired injection conditions are reached. 

As mentioned in the section TRAJECTORY DESCRIPTION, the trajectories consist 
of a short vertical rise followed by a rapid pitchover phase determining the amount of tra- 
jectory lofting during the near zero angle-of-attack atmospheric phase. This near zero 
angle-of-attack phase lasts until the booster engines are jettisoned. Then the steering is 
determined by the variational solution. 

The first mission is a sun-synchronous circular orbit at 1110 kilometers altitude. 

An inclination of 99. 9° is required to take advantage of precession due to oblateness to 
maintain sun synchronization. A launch azimuth of about 193° would be required for a 
planar trajectory to the required orbit. Atlas-Centaur trajectories from ETR at that 
launch azimuth present severe, if not insurmountable, range safety problems because 
scheduled hardware drops (booster engines, insulation panels, payload fairing, Atlas 
stage) as well as the HP trace endanger land areas. Hence, a launch azimuth of 108° was 
chosen for the demonstration trajectory to avoid these hardware impact problems. The 
thrust attitude is constrained to be parallel to the launch azimuth plane until 5 seconds 
after Centaur ignition (268 sec after lift-off) after which the attitude is unconstrained. 

The rate of change of attitude is constrained to be less than or equal to 1. 6° per second. 
The free parameter during the atmospheric portion, the pitchover, is optimized. The 
thrust attitude from booster jettison to injection is optimized using the analysis and tech- 
niques described in this report. The HP trace for this trajectory is presented in figure 1. 
A spherical earth model with no drag is used to calculate the IIP trace. The injected 
payload weight is 800 kilograms. A trajectory with the same ground rules, but with no 
attitude rate constraint, achieved an injected payload weight of 820 kilograms. Thus, 
the payload loss due to the attitude rate limit is 20 kilograms. The two HP traces are 
virtually identical. The thrust attitude in pitch (ip) and yaw (j3) (see fig. 2) are presented 
in figures 3 and 4, respectively. The thrust attitude rate is dis continuously increased 
to the limiting rate of 1.6° per second 0.07 second before the thrust attitude is freed 
from the planar constraint and remains at that level for the following 34 seconds. Fig- 
ures 3 and 4 show that during the planar portion of the trajectory the thrust vector is 
directed such that altitude is gained at the expense of horizontal velocity. It is the hori- 
zontal velocity which must be turned to satisfy the final inclination. After the planar 
constraint is removed, the thrust attitude rapidly approaches the horizontal to accumu- 
late the required circular velocity. The total rate of change of thrust angle, altitude, 
inertial flight path angle, and velocity as functions of time are presented in figures 5, 6, 

7, and 8, respectively. 

Launch trajectories from ETR to sun-synchronous orbits unavoidably present range 
safety problems because the instantaneous impact point traces will always cross popu- 
lated land areas in the Caribbean and/or South America. It is impossible to avoid over- 
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fly for any practical trajectory regardless of the method used to alter the HP trace. The 
trajectory chosen in this report does not necessarily reflect the best choice from a range 
safety standpoint. It is used only to demonstrate an application of the analysis to a mis- 
sion where significant attitude and attitude rate constraints are present. 

The second demonstration case is characteristic of another range safety problem 
whereby the HP trace can be moved to avoid all land masses. For certain combinations 
of declination and right ascension of the outgoing geocentric asymptote, planar inter- 
planetary trajectories possess impact point traces which cross South America, in parti- 
cular Venezeula and Brazil. It is possible in some of these cases to move the IIP trace 
north to avoid all land masses. This is demonstrated in the trace presented in figure 9. 
The declination and right ascension of the outgoing asymptote are -45° and 190. 5° east of 
the launch site meridian, respectively. The injection energy is 4.0215 square kilometers 
per square second with a perigee altitude of 167 kilometers. The planar trajectory would 
require about a 200° launch azimuth. A planar trajectory at this launch azimuth would 
be unacceptable because the hardware impact points and IIP traces impinge land. How- 
ever, land masses are completely avoided by choosing a launch azimuth of 108° and by 
constraining the thrust attitude to be parallel to the launch azimuth plane until 150 seconds 
after Centaur ignition. The rate of change of thrust attitude is again limited to 1.6° per 
second. The pitchover phase is optimum. The thrust attitude rate is increased discon- 
tinuously to the limiting value 0.4 second before the planar constraint is removed and 
continues at the maximum rate for the following 17 seconds. The injected payload weight 
for this case is about 227 kilograms. A similar trajectory without the attitude rate limit 
constraint achieved about a 2 kilogram increase in payload weight. Figures 10 to 15 
show the thrust angles in pitch and yaw, the total rate of change of the thrust angle, the 
inertial flight path angle, and the velocity and altitude, respectively, as functions of time. 


CONCLUDING REMARKS 

The problem of obtaining the required terminal conditions for orbital and planetary 
missions while satisfying typical range safety constraints has been considered. To 
simplify the mathematical problem, a planar attitude constraint was imposed for the 
initial portion of the flight in order to avoid overflight of certain land areas. Later in 
flight, the attitude constraint is removed in order to allow the required terminal condi- 
tions to be achieved. An attitude rate limit is also imposed in order to prevent discon- 
tinuous attitude changes. 

The mathematical problem is formulated by using the calculus of variations. Euler- 
Lagrange equations are obtained which determine the optimum thrust direction profile 
subject to the problem constraints. A simplified procedure is also presented for obtain- 
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ing the auxiliary variational boundary conditions required for free end conditions. 

To illustrate the use of the equations, results are presented for two typical missions 
which require nonplanar trajectories of the type discussed. These missions are a sun- 
synchronous orbit and a typical interplanetary probe. An Atlas -Centaur launch vehicle 
is assumed, launched from the Eastern Test Range. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 3, 1969, 

128-31-30-02-22. 
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APPENDIX A 


SYMBOLS 


a 

o 

thrust acceleration, m/sec 

s 

sign (f x a) • n 

b i 

constraint equations 

t 

time, sec 

b j 

constraint equations 

V 

velocity vector, m/ sec 

c 

first integral of Euler -Lagrange 

x i 

problem variables 


equations 


constant vector 


y 

d() 

differential of ( ) 


unit vector pointing at North Pole 

z 

E 

Weierstrass excess function 

a 

arbitrary scalar 

e 

eccentricity 

0 

yaw attitude, deg 

F 

defined in eq. (lb) 

7 

Lagrange multiplier, sec 

f 

thrust direction vector 

A() 

jump discontinuity in variable 

G 

Earth gravitational constant. 

5() 

denotes linearized variable 


3/ 2 

m /sec 





e 

arbitrary vector 

G 

gravity acceleration vector, 



m/ sec 2 

V 

true anomaly, rad 

g 

function of initial and final con- 

9 o 

commanded attitude discontinuity, 


ditions to be minimized, sec 


rad 

H 

2 2 

energy, m / sec 

e 

unit tangential vector 

h 

2 

angular momentum, m /sec 

X 

2 

Lagrange multipliers, sec /m 

I 

orbital inclination, rad 

M 

Lagrange multipliers, sec/m 

J 

functional to be minimized by 

p 

Lagrange multiplier, sec 


variational methods, sec 

CT mag 

A 

magnitude of f x a at t^, sec 

m 

mass, kg 

a 

Lagrange multipliers, sec 

n 

unit normal to attitude con- 

<P 

angle from injection to outgoing 


straint plane 


asymptote, rad 

P 

semilatus rectum, m 

* 

pitch attitude, rad 


generalized state variables 

w m 

maximum rotation rate, rad/ sec 

r 

radius, m 

CO 

rotation vector, rad/sec 

S 

constant vector in eq. (B6) 
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Subscripts: 

1 

start of regular planar arc 

A 

asymptote 

2 

start of unconstrained phase 

av 

average value 

3 

start of singular unconstrained 

d 

desired 


arc 

f 

final 

oo 

hyperbolic excess 

h 

normal 

Superscripts: 

I 

injection 

• 

derivative with respect to time 

i 

problem variables 

* 

finite, but allowable variation of 

P 

perigee 


the variable 

r 

radial 

*A> 

unit vector 

t 

fixed time 

— 

vector 

9 

tangential 

+ 

evaluated after switching point 

0 

initial 

- 

evaluated before switching point 
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APPENDIX B 


EQUATIONS FOR PROBLEM SOLUTION 


The variational portion of the trajectory is divided into four phases, as discussed 
earlier. These phases, along with the definition of switching times, are presented in 
table I. 

The following equations must be integrated throughout the variational portion of the 
trajectory (eqs. (2a), (2b), (3a), and (3b)). 


v + GmT-af=0 




r - v = 0 

X + jJL = 0 


> 


£-Gmi + 3 ^(X • 7)7 = 0 

3 5 

r r 


(Bl) 


These equations are linearized, and the linearized equations are also integrated through- 
out the trajectory in order to obtain the partial derivatives required for the Newton- 
Raphson iteration. The linearized equations obtained from equations (Bl) are 


6 v 



3 < ^ m (r • 6r)r - a6? = 0 

5 

r 




6r - 6v = 0 
6 A + 5 M = 0 


> 


67T - 2“6X + Rr • 6r)X + (6A • r)r + (X • 6r)r + (X • r)6r"| - 15^ (X ■ r)(? • 5r)r =15 


(B2) 


where 6 is used to indicate a linearized variable. As noted in the Numerical Iteration 
section, fxa is not required until phase 2. It will be shown later that a closed form 
solution is available for f x cr during phase 3. Also, f Xa is identically zero during 
phase 4. The differential equation for f must be integrated during phases 2 and 3, while 
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the singular solution is used during phases 1 and 4. The additional equations to be inte- 
grated are (from eqs. (2c) and (3c)) 


Phase 2 


— (f x a) = aX x f - a) x (a x f ) 
dt 


"\ 


— fs (f x a)l = a6X x f + aX x 6f - o> x 5(0 x f) 
dt L J 




A A 

f = ( 1 ) Xf 


6f = co x 5f 




where u> = sign [<f x a) ■ njn (eq. (9)). 


Phase 3 


f = O) xf 


5f = at x 5f + 6w x f 


(B3) 


(B4) 


where (from eq. (8)) 


co 

m |f X a 


and 


5(f X g) _ (f X cr)(f X q) • 5(f X g) 


During phase 3, f Xa can be calculated in closed form by using the fact that 

AXv + fiXr-fXcf: constant (B5) 

Equation (B5) can be verified by differentiation and using equations (2), (3), and (10). 
Therefore, f x a can be calculated from 
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fXCT = XXv + juXr-S 


(B6) 


where 


S = (XXv + /uXr-fXa). 

*2 

Taking the variations of equations (B6) and (B7) results in 

6 (f X a) = (6 X X v + X X 5 v + 8 Jl Xr+]IX5r- 6S) 


and 


6S = |6XXv + XX6v + 5juXr + /i X6r-6 



Calculation of Thrust Direction 

During phase 1, the thrust direction is calculated from equation (14): 


£ _ X - (X. • n)n 

Vl - (* • n) 2 


The variation of equation (B10) is 

g£ _ 5X - (5X • n)n + |]x - (X • n)n](X • n)(5X • n) 

Vi - ft • S ) 2 [i . (X . S)f /2 

where 


6X = — - — (6X • X) 
X x 3 


During phase 4, f is calculated from equation (12): 

f = X 


and 


(B7) 


(B8) 


(B9) 


(BIO) 


(Bll) 
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Sf = — - -L (6X • X) 

X x 3 

The thrust direction is integrated during phases 2 and 3, as indicated earlier. 


Initial Conditions 

The choice of initial conditions is discussed in appendix C, and the variable initial 

conditions are i//, i|/, j3, /3, and X. In addition, the switching times t^ and tg are used 

as variable initial conditions along with the value of cr at t, . Equations must be 

ma &_ _ 1 ^ _ 

derived for the linearized changes in the initial values of X, p, and f x a, obtained as a 
result of linearized changes in the variable initial conditions. 

The initial variations in X and p are obtained by linearizing equations (C5b) and 
(C5c): 


5X = X5% 

5p = -X5X - 6 XX - X5X 

The initial value of X is always equal to unity, and the initial values of r , §, and fi are 
fixed because the initial position and velocity are specified. Also, the initial value of ? 
is fixed because r is a function only of position and velocity and not thrust direction. 
Therefore, the variations of equations (C5a) and (C6) are 

5X = 5X r r + 5X e 6 + 6X h h 


and 


6X — 5Xj,r + 5X r r + &Xq6 + 5Xgd + Xq58 + 5X^h + 5Xj^h + X^5h 


Equations (Cl), (C2), and (C3) are linearized to give 


6X r = -sin p sin ^6p + cos p cos tybip 


5 Xq = -sin p cos 5/3 - cos p sin i//5 xf/ 


6X h = cos pep 


26 



5X r = -sin j3 sin \p6(3 + cos /3 cos \pb\p - (/3 cos /3 sin xp + xp sin 0 cos i//)60 

-(0 sin 0 cos xp + xp cos 0 sin 

6X 0 = -sin 0 cos ^60 - cos 0 sin i^5)// - (0 cos 0 cos xp - xp sin 0 sin j0)60 

+(0 sin 0 sin xp - xp cos 0 cos 
6X h = cos 060 - 0 sin 06/3 


60 = (8f • h)h 

h 


6h = - H. ( 6 ! • fi)fi 

h 


where 6f is calculated by using equation (Bll). 

Some of the linearized variables have jump discontinuities at a definite time due to 
the allowable variation in this time. These jumps occur when the derivative of the vari- 
able is discontinuous at this time and are calculated from equations of the form 

(5x)* = (6x)" + (x £ - xj^St (B12) 


where 6x^ is the linearized change at fixed time and the minus and plus superscripts 
denote evaluation of the variable before and after t, respectively. 

The derivative of f is discontinuous at tj. The derivative after t^ is given by 
equation (9). 



co x f = sign 


[jf X a) • njn X f 


and (f). 

H 


is evaluated by differentiating equation (BIO): 



X - (X • n)n + [x - (X • n)n](X • n)(X • n) 

» /j -<* ■ s > 2 [i - <x • i) z f /2 


where 
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£ _ X _ x(x • X) _ _JU X(X • /i) 

X x 3 XX 

Therefore, by using equation (B12), 

(81)* = (Sf)^ + (f* - 6tj 

The derivative of f x cr is also discontinuous at tj, but the situation is somewhat 
different than for f, since f x a is initialized at tj. When f x <j is initialized at tj, 
the value is (from eq. (27)) 

The value of f x a at a later time (tj + 5tj) is obtained by using equations (2c), (3c), 
and (9). 

< f x + 8t 1 ) = <W (f x ft + 1 <f * ^>«i 


= cr mag (f X n) + a(X X fJStj - co X (^ X f)6t x 
= <J mag (f x n) + a(X X j^etj + < 7 mag w m sn x (? x fi)5t 1 (B12) 


where 


s = sign(? x a) • n 

When f x a is initialized at a later time (tj + 5tj), the value is 

^ x a) (t x + 6t x ) = a mag ( * X n) (t x + 5t x ) 


where , by using equations (2c) and (9) 


(f x n) 


(tj + 6t x ) 


.A 

(f x n). 
r l 

(f x n), 
t l 


+ (co x f) x n5tj 
+ ct» m s(n x f ) x n5tj 
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Then f x a becomes 

(? X ^(t x + Stj) = CT mag( f X “> + cr mag W m 8 < S X X " 6t l < B13 > 

The jump discontinuity in f X a is given by the difference between equations (B13) and 
(B12): 

A(f Xu) = -a (X Xfje^ 

and the variation of f x <j at tj is given by 

8(f xj) = 6o mag (f X n) + u mag (5f X n) -a(X X fJStj 

During the iteration, the value of f is allowed to be discontinuous at tg. This 
results in a jump discontinuity in 5v. The total 5v is given by 

(5v)^ = (6v)" + a(f - X)“ 5t„ 

Final Conditions 

The variational final conditions for the sun-synchronous mission are (see eq. (27)). 

(f x a). =0 

r 3 

, a* ^ . 

(f * X) t =1 
*3 

(X x v + ju x r) f • z = 0 

t f 

|(X X v + ju, x r) • h)jj. = 0 

For the interplanetary mission, the variational final conditions are 

(f X <t), =0 

t 3 
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(f • X) t =1 
x 3 

(X • G + M ' v). =0 

t f 

(X x v + ju x r) • =0 

h 

In addition, there are three final conditions on state variables for both of the missions 
(energy is used for termination of the trajectory). The final state variable end conditions 
for the planetary mission are 

(1) Perigee radius 

(2) Declination of outgoing asymptote 

(3) Right ascension of outgoing asymptote 

Conditions 2 and 3 are equivalent to specified v w . For the sun-synchronous mission, the 
state variable end conditions are 

(1) Radius 

(2) Radius rate 

(3) Inclination 

In order for the iteration to converge on these final conditions, the linearized changes 
in final conditions are required. The changes in the variational end conditions at tg are 
obtained by using equations (2c), (3a), (10), and (B12). 


[5(fxaj t + = [5(f X a)]~ -a[(x x f j| t " 6t 3 

3 3 3 



where u> is evaluated by using equation (8) . The changes in the variational final condi- 
tions at tj, are 


6(X • G + ju • v) t = 




i -n i Rii . XT _L I I . F\ir\ 


Sr + 3 (X • r)(r 
r 5 


Sr) + Sju • v + ju • Sv 


v 

- 


(B14) 


i|(xxv+ pcxr) • x] t = |§Xxv + XX6v + juxfir + fiptxr)]^ • x^ + (XXv + MXr) t ^ 


6x 
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The 6r , 6v, 6A, and 5~jl at t f in the preceding equations have been modified to be 
variations at constant energy rather than constant time. The following equations are used 
to modify these variables : 


6t, 


SH 

H 



Sv r • 5r 

r 3 _ 
aA • v 


-It, 


6r = (5r) t + v^6t f 
5v = (6v) t + |aA - r j 6t f 


SA — (SA)^, — 5tj 


5/i = (6ju.) tf + M tf 6t f 


where H is the injection energy 

H = ^ ' Z - G m 
2 r 

The 6x in the last of equations (B14) is zero if H is "v^ or z. If x corresponds to h, 
then 


6x = 6h = (5r Xv + rX 5v) f 

i 

In order to determine the linearized changes in final state conditions, these final 
conditions must first be expressed in terms of r and v. The perigee radius is expres- 
sed as follows : 


r =-P- 
p 1 +e 
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where 


h^ _ (r x y) • (r x y) 
Gm Gm 


and 



The direction of the outgoing asymptote is given by 

= r cos <p + 8 sin <p 

(see sketch b), where <P = ?7 A - and the true anomalies are given by (ref. 7) 



Vco 

■^Outgoing 

asymptote 


V A = cos 


X-;) 



77j = sin 


The radius rate and inclination (for the sun-synchronous mission) may be calculated 
from 
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r = r • v 


I = cos' 


‘(¥) 


The linearized changes in the state variable end conditions are then given by 


6r„ 


= 2 l h ’ r (l - -i!L\ 

h 2 P \ 


5v oo = 5r cos ^ + & sin (p + (6 cos ^ - r sin <p)5<p 


where 


= 6r _ (r • 5r) - 


6r = — 


r r 3 


6 fi - Oh (h • 5h)h 


h h 3 


86 = 6h x r + h x 6r 


8V A = 


8(p = 8t] A - 5 t?j 
2Hh • 6h 


(Gm) 2 e 2 Ve 2 - 1 


1 


O^^tan^l b ' g . h +. 8r ' v + r • 6v 2Hh • 6h r • dr' 


Gmp 


r • v 


(Gm) 2 e 2 r 2 


The linearized changes in r and I are given by 


6r = 6r • v + r • 5v 


61 = - 


5h • z 
sin I 
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APPENDIX C 


CHOICE OF INITIAL CONDITIONS 


Solution of the two point boundary value problem requires initial guesses for eight 
variables, six of which are Lagrange multipliers. These variables are listed in equa- 
tion (27). In order to eliminate the difficulty in guessing at values of these multipliers, 
equations are derived which express four of the multipliers in terms of pitch and yaw 
attitude and rate. In addition, equations are presented for deriving reasonable first 
guesses of the other multipliers and switching times. 

The pitch attitude xp and yaw attitude /3 are defined in figure 2. If the attitude were 
not constrained to be planar, the vehicle attitude would be along the X vector. There- 
fore, the guesses of xp, xp, (3, and /3 correspond to the optimum values for the uncon- 
strained attitude case. The coordinate axes in figure 2 correspond to a unit radial, tan- 
gential, normal coordinate system defined by 


- ^ 
r = — 


h = 


r x y 
r x v 




$ = fi x r 


(Cl) 


The following equations are easily obtained from the figure: 


X r = X • r = cos /3 sin xp 


\q = X • 6 = cos jS cos xp f 


X. = X • h = sin /3 


(C2) 


Equations (C2) are differentiated to give 

X r = -/§ sin /3 sin xp + xp cos /3 cos i p'' 

sin /3 cos xp - xp cos /3 sin xp > 

X h = 0 cos p 


(C3) 
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The coordinate system defined in equations (Cl) is not inertial. The derivatives of 
equations (Cl) are given by 






A? V 


h 


h = - H x. 


J 


(C4) 


where 


h = |r x v 


The values of X and p. can now be calculated from 


A -A A A 

X = X r r + Xg0 + X^h 


X = XX 


/! = -X = -XX - XX 
where X is calculated from equations (C3) and (C4) by using 


X = X r ? + \q6 + X^h + X r ? + \q6 + X^fi 


(C 5) 


(C6) 


It remains to supply values for X and X in order to completely determine X and /i . 
Recalling that the Euler Lagrange equations are homogeneous in the multipliers, X can 
be set equal to unity without loss of generality. In order to estimate X, it is first noted 
that if there is no planar attitude constraint, the vector 


XXv + jLiXr + aXf 


(C7) 


is a constant of the motion. This can easily be verified by differentiating and using 
equations (2) and (3). Both of the missions discussed herein contain required final con- 
ditions of the form 
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(C8) 


(xxv + juxr)-y = 0 


(see eq. (20) and (22)) where y is some fixed vector (z or v TO ). Also, if the initial and 
final attitudes are not specified, a x ? = 0 at the beginning and end of the variational 
problem. Therefore, under these restricted conditions, equation (C8) can be satisfied at 
the beginning of the trajectory, thus reducing the iteration size. Equation (C8) is satisfied 
by the proper choice of the initial value of X. Combining equations (C5) and (C8) results 
in 


X = X (^ x v - A x r) • y 
(X x 7 ) • y 


(C9) 


Because of the attitude constraint during part of the trajectory, the vector in equa- 
tion (C7) is not constant throughout the trajectory. Therefore, equation (C9) is not exact 
and the iteration size cannot be reduced. However, equation (C9) still supplies a reason- 
able first guess for X. 

A reasonable guess for tg may be obtained by estimating the angular discontinuity 
at tg for the nonrate- limited case and dividing by the maximum vehicle turning rate: 


cos -1 (X • f) tf 


t 3 =t 2 + 


o> 


m 


The value of tj is more difficult to estimate, but (tg - tj) should generally be an order 
of magnitude smaller than (tg - tg). 

An additional multiplier remains to be estimated, which is the value of CT ma g in the 
equation 


< f x a >tj = "mag 9 * X 

The estimating procedure for a ma g starts by noting that 

|fx " 1 t|=(‘w) ti 


and | f x a 


must be zero at tg 


Also, 


— (f x a) = -af x x 
dt 
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between t g and tg. Since the time interval Oq, tg) is much shorter than (tg, tg), v mSL g 
may be estimated as 


a 


mag 



dt 


(CIO) 


During the interval (tg , tg) , the motion of £ is essentially toward X, while X remains 
nearly constant. Therefore, the direction of fxx is essentially fixed, so that 


— If x cr I « -alf x xl 
dt 


and 



t Q and i 
becomes 


(Cll) 


where 


The value of X in this integral may be assumed to be unity, since X = 1 at 
is very small. If a is assumed constant (equal to its average value), cr 


mag 


mag 


/ t 3 -t 2 


a sin (0 - a) t) dt 

' o m 7 


av 


00 


(l - cos e» o ) 


m 


cos e Q = (X • f) t2 


and a ay is the average acceleration between tg and tg. 
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TABLE I. - TRAJECTORY PHASES AND SWITCHING TIMES 


Event time 

Description 

Sun-synchronous 

trajectory 

Interplanetary 

trajectory 



Time of occurrence, sec 

Liftoff 

Vertical rise 

0 

0 


Pitchover and near-zero 
angle of attack booster 
trajectory 

20 

20 

*0 

Start of singular planar 
arc (phase 1) 

156.00 

155.27 

h 

Start of regular planar 
arc (phase 2) 

267.73 

405.44 

t 2 

Start of regular uncon- 
strained arc (phase 3) 

267.80 

405.86 


Start of singular uncon- 
strained arc (phase 4) 

301.76 

422.45 


Powered flight termin- 
ation 

706.54 

710.06 











1 2 3 4 5 6 7 

Time of flight, sec 

Figure 11. - Vehicle yaw attitude as function of time. Interplanetary trajectory. 



1 2 3 4 5 6 7 

Time of flight, sec 

Figure 12. - Total vehicle turning rate as function of time. Interplanetary trajectory. 
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